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ABSTRACT 

The forkhead transcription factor FoxOI is a critical 
regulator of hepatic glucose and lipid metabolism, 
and dysregulation of FoxOI function has been 
implicated in diabetes and insulin resistance. We 
globally identified FoxOI occupancy in mouse 
hepatic chromatin on a genome-wide level by 
chromatin immunoprecipitation coupled with high- 
throughput DNA sequencing (ChlP-seq). To estab- 
lish the specific functional significance of FoxOI 
against other FoxO proteins, ChlP-seq was per- 
formed with chromatin from liver-specific FoxOI 
knockout and wild-type mice. Here we identified 
401 genome-wide FoxOI -binding locations. Motif 
search reveals a sequence element, 5 GTAAACA 
3 , consistent with a previously known FoxOI - 
binding site. Gene set enrichment analysis shows 
that the data from FoxOI ChlP-seq are highly 
correlated with the global expression profiling of 
genes regulated by FoxOI, demonstrating the func- 
tional relevance of our FoxOI ChlP-seq study. 
Interestingly, gene ontology analysis reveals the 
functional significance of FoxOI in retinoid meta- 
bolic processes. We show here that FoxOI directly 
binds to the genomic sites for the genes in retinoid 
metabolism. Notably, deletion of FoxOI caused a 
significantly reduced induction of Pck1 and Pdk4 in 
response to retinoids. As Pck1 and Pdk4 are down- 
stream targets of retinoid signaling, these results 
suggest that FoxOI plays a potential role in linking 
retinoid metabolism to hepatic gluconeogenesis. 



INTRODUCTION 

The liver plays a central role in maintaining tight control 
of glucose homeostasis. Under feeding conditions, there is 
increased hepatic glucose uptake and glycogen storage, 
whereas during fasting, the liver produces glucose by 
gluconeogenesis and glycogenolysis (1). Glucose homeo- 
stasis during this transition period is orchestrated by 
the two counteracting pancreatic hormones, insulin and 
glucagon. During feeding, insulin activates phospho- 
inositide 3-kinase (PI3K)/Akt (also known as PKB) sig- 
naling, which induces phosphorylation of forkhead 
transcription factor FoxOI, leading to the export of 
FoxOI to the cytoplasm from the nucleus (2). As a 
result, insulin inhibits expression of FoxOI target genes, 
including those encoding phosphoenolpyruvate carboxy- 
kinase (PCK1) and glucose-6-phosphatase (G6PC), the 
rate-limiting enzymes in hepatic gluconeogenesis (3,4). 

FoxOI is a member of the evolutionarily conserved 
mammalian forkhead box O class (FoxO) family that also 
includes Fox03, Fox04 and Fox06. FoxO proteins bind to 
their consensus DNA sequence, (T/C)(G/A)AAACAA, as 
a monomer, and regulate expression of genes involved in a 
variety of biological processes, including cell cycle arrest, 
apoptosis, detoxification, energy metabolism and longevity 

(5) . Studies using genetically modified mouse models have 
well documented the physiological role of FoxOI in vivo in 
insulin action in suppressing hepatic glucose production 

(6) . Additionally, recent studies have shown that the in- 
volvement of Notch (7) and Wnt (8) signaling in regulating 
hepatic glucose production is through modulating FoxOI 
activity. There is also an emerging theme suggesting a role 
of FoxOI in hepatic lipid homeostasis. Consistent with 
its critical role in hepatic metabolism, aberrant function 
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of FoxOl is associated with deleterious characteristics of 
type 2 diabetes, including hyperglycemia and hypertrigly- 
ceridemia (9). 

Given the critical function of FoxOl in hepatic metab- 
olism, it is necessary to completely understand how 
FoxOl functions in controlling the diverse metabolic 
processes. To this end, we have performed genome-wide 
analysis of FoxOl -binding occupancy in hepatic chroma- 
tin by chromatin immunoprecipitation coupled with 
high-throughput DNA sequencing (ChlP-seq) (10). Here 
we have identified 401 FoxOl -binding locations in mouse 
hepatic chromatin. Gene ontology (GO) analysis reveals a 
significant role of FoxOl in metabolic processes for 
carboxylic acids, fatty acids and steroids. Interestingly, 
we also observed a cluster of genes involved in retinoid 
(vitamin A) metabolic processes. 

Retinoids (retinol derivatives) play critical roles in de- 
velopment, immune function, vision, tissue modeling and 
metabolism. The liver is the major organ for retinoid 
uptake and storage in the body (11,12). Biologically 
active retinoid metabolites include retinoic acids that are 
generated by a two-step process from retinol, catalyzed by 
retinol dehydrogenases and retinaldehyde dehydrogen- 
ases. All-trans retinoic acid (atRA) is a transcriptionally 
active metabolite that functions as a ligand for its specific 
receptors, retinoic acid receptors (RARs), and directly 
regulates its target genes with its heterodimeric partners, 
retinoid X receptors. Cellular retinol-binding protein 1 
encoded by retinol-binding protein 1 (Rbpl) serves as a 
carrier protein for retinol and retinal. 

It has been shown that retinoids play an essential role 
in hepatic glucose metabolism by regulating expression of 
gluconeogenic gene Pckl (13-16). Also, pyruvate dehydro- 
genase kinase isozyme 4 (Pdk4) expression levels are 
regulated by vitamin A in vivo (17,18). Therefore, we were 
interested in further studying the role of FoxOl in retinoid 
metabolic processes and its downstream effects on hepatic 
gluconeogenesis. We show here that FoxOl is required for 
full induction of Pckl and Pdk4 in response to retinoids, 
suggesting a previously unappreciated role of FoxOl in 
retinoid signaling, regulating hepatic gluconeogenesis. 

MATERIALS AND METHODS 

Animals 

All animal care and experimental procedures were per- 
formed according to the guidelines of the Institutional 
Animal Care and Use Committee. Liver-specific FoxOl 
(1-FoxOl) knockout (KO) mice were produced by 
crossing Fox01 flox/flox mice expressing CRE recombinase 
under the control of the ocl -antitrypsin promoter, as 
described in a previous study (19). Littermates FoxOl flox/ 
flox were used as wild- type (WT) controls. Mice were main- 
tained on a 12-h light/dark cycle with free access of food and 
water. Eight- to 9-week-old male mice were either fed ad 
libitum or fasted for 24 h. Livers were removed and pro- 
cessed immediately for chromatin isolation, as described 
previously. The remaining 10% of liver was frozen in 
liquid nitrogen and stored at — 80°C until RNA and 
protein were isolated. 



ChIP and ChIP sequencing 

Mouse liver chromatin was prepared, as described previ- 
ously (19). An aliquot of chromatin was pre-cleared with 
protein A and immunoprecipitated with an anti-FoxOl 
(Cell Signaling Technology). ChlP-enriched DNA was 
isolated and further purified by QIAquick PCR purifica- 
tion kit (Qiagen). After verification of FoxOl enrichment 
to known FoxOl target promoters by quantitative real- 
time polymerase chain reaction (qPCR) using gene-specific 
primers, 20 ng of ChlP-enriched DNA was sent to Ambry 
Genetics (Aliso Viejo, CA), where it was processed for 
high-throughput DNA sequencing. Briefly, the DNA frag- 
ments were blunted and ligated to sequencing adapters. 
Samples were size-fractionated on a 2% agarose gel, and 
DNA fragments 200 ± 25 bp in length were isolated. After 
PCR amplification, the resulting DNA library was 
sequenced on the Solexa/Illumina Genome Analyzer. 
The sequence reads of 36 bases were aligned to the mm9 
mouse genome using ELAND, allowing one mismatch. 
The aligned reads were used as input to identify peaks 
as well as visualization in the UCSC Genome Browser. 
To compare the genomic locations of our FoxOl ChlP- 
seq peaks with those of FoxA2 ChlP-seq peaks, the 
FoxOl ChlP-seq bed files aligned to the mm9 mouse ref- 
erence genome were re-mapped onto the mm8 mouse ref- 
erence genome using UCSC liftOver (http://genome.ucsc. 
edu/cgi-bin/hgLiftOver), as the sequence reads of FoxA2 
ChlP-seq were aligned onto the mouse reference genome 
version mm8. The FoxA2 ChlP-seq data were obtained 
from http://www.bcgsc.ca/data/chipseq/chip-seq-data 
(20). 

Peak finding, density map plot and motif analysis 

Peak finding was performed to identify significantly 
enriched reads in peaks in WT chromatin ChlPed with 
an anti-FoxOl relative to those in WT immunoglobulin 
G (IgG) or in 1-FoxOl KO chromatin ChlPed with an 
anti-FoxOl by using SEAN, custom peak finding 
software developed by the University of Connecticut 
Bioinformatics and Bio-Computing Institute 
(Supplementary Materials and Methods SI) with a 
default P- value cut-off criterion of 1 x 10 -5 to generate 
a high-confidence peak set. Density map plots were 
produced from each read in peaks 5kb upstream and 
downstream from the center of the peak. Each read was 
displayed horizontally with a green spot. Density maps for 
random peaks were produced from reads in the peaks of 
randomly selected lOkb of genomic regions. 

To search for sequence motifs representing FoxOl - 
binding sites, the underlying DNA sequences of the 
identified peaks were obtained by Galaxy (main.g2.bx. 
psu.edu) and fed into motif-finding program MEME 
4.7.0 with default parameters (21). Motifs were repre- 
sented as position-dependent letter-probability matrices, 
which describe the probability of each possible letter at 
each position in the pattern. To identify enriched 
binding sites for co-regulatory proteins in FoxOl peaks, 
the regions spanning 250 bp upstream and downstream 
from the FoxOl motifs were scanned by HOMER motif 
analysis. 



Nucleic Acids Research, 2012, Vol 40, No. 22 11501 



Peak distribution, GO analysis and gene set enrichment 
analysis 

The peak distribution was computed based on the UCSC 
refFlat file database. The distance from the peak to the 
nearest transcription start site (TSS) was analysed and 
plotted. The random peaks were placed on the genome 
at random locations and determined the distance to 
TSS, and a background plot was produced. The TSSs 
were referenced from the UCSC refFlat file. GO analysis 
was performed with nearest genes of FoxOl -binding 
peaks as input to the National Institutes of Health 
Database for Annotation, Visualization and Integrated 
Discovery GO (22). Gene set enrichment analysis 
(GSEA) was carried out by comparing a data set of 
FoxOl ChlP-seq genes with that from genome-wide 
mRNA expression profiling (23). Enrichment score was 
obtained with a P-value by running sum of the ChlP-seq 
genes and differentially expressed rank-ordered genes in 
mRNA profiling. 

RNA isolation, qPCR and microarray data analysis 

Total RNA was isolated from mouse livers or primary 
mouse hepatocytes using TRIzol (Invitrogen) according 
to the manufacturer's instructions. cDNA was synthesized 
using the iScript cDNA synthesis kit (Bio-Rad). qPCR 
was performed with the iQ5 real-time PCR detection 
system (Bio-Rad), and the results were analysed by using 
the AA-Ct method supplied by the manufacturer, as 
described previously (19). mRNA levels were normalized 
to the level of mRNA for the housekeeping ribosomal 
protein L32. Gene expression profiling was performed 
using the Mouse Gene 1.0 OST (Affymetrix, Santa 
Clara, CA) by hybridizing RNA from fed or fasted 
livers of WT and fasted 1-FoxOl KO livers in triplicate. 
cRNA generation, labeling and hybridization were per- 
formed at the DNA Microarray Facility of the 
University of California, Irvine. Expression values were 
analysed using Limma (Linear Models for Microarray 
Data) (24), SAM (25) and Bayesian-based 
Semi-Parametric Hierarchical Model (26) for statistical 
significance. 

Primary mouse hepatocytes, adenovirus infection and 
retinoid treatment 

Primary mouse hepatocytes were isolated from 8 to 10 
weeks old C57BL/6J male mice (The Jackson 
Laboratory), as described previously (19). Briefly, hepato- 
cytes were obtained by two-step collagenase perfusion 
through the cannulated portal vein, followed by Percoll 
gradient selection. Isolated cells were resuspended in 
medium 199 supplemented with 10% FBS, 10 nM dexa- 
methasone, 2mM glutamine, 100 U/ml penicillin, 
lOOU/ml streptomycin and 1 x ITS, and plated in 
collagen-coated 12-well dishes at a density of 2.5 x 10 5 / 
well. After 4-h attachment, the medium was replaced 
with serum-free medium 199 supplemented with lOnM 
dexamethasone, 2mM glutamine, 100 U/ml penicillin 
and 100 U/ml streptomycin. Sixteen hours after isolation, 
hepatocytes were infected with a recombinant adenovirus 



expressing un-speciflc RNAi or FoxOl shRNA (gift from 
Dr. D. Accili, Columbia University) at a multiplicity of 
infection of 25 and incubated for 24 h in medium 199 sup- 
plemented with 2% serum, followed by treatment with 
retinol (Sigma) or atRA (Sigma). Six hours later, cells 
were harvested for RNA analysis. 

RESULTS 

We prepared chromatin from mouse livers after 24-h 
fasting because FoxOl is predominantly localized in the 
nucleus during fasting, promoting the expression of its 
target genes. Before carrying out ChlP-seq analysis, we 
first evaluated the quality of an antibody against FoxOl 
in liver- specific FoxOl (1-FoxOl) KO mice. Immunoblot 
analysis showed (Supplementary Figure SI A) that hepatic 
FoxOl protein was barely detectable in 1-FoxOl KO, 
compared with WT littermates. The protein levels were 
proportionally reflected with mRNA levels of FoxOl 
that were significantly decreased more than 90% in 
1-FoxOl KO (Supplementary Figure SIB). We used this 
antibody for gene-specific ChIP with chromatin from WT 
and 1-FoxOl KO and analysed the enrichment of FoxOl 
to known FoxOl -binding sites to verify the ability of the 
antibody to immunoprecipitate FoxOl. We observed sig- 
nificant enrichment of FoxOl at the promoter region of 
the Pckl gene, which has been previously shown occupied 
by FoxOl (3), compared with IgG in WT, whereas a 
genomic fragment of Yyl locus used as a negative 
control was not enriched (Supplementary Figure S2). 
The significant binding of FoxOl in WT was further 
verified by ChIP with chromatin from 1-FoxOl KO 
where there was no specific FoxOl binding to the Pckl 
promoter regions observed (Supplementary Figure S2). 
Consequently, we used the antibody for ChlP-seq for a 
genome-wide analysis of FoxOl binding. Because there 
are several members of FoxO proteins that share a 
DNA-binding domain, in addition to an IgG control, we 
prepared an additional negative control for the specificity 
of FoxOl against other FoxO proteins. Hepatic chroma- 
tin from 1-FoxOl KO was immunoprecipitated with a 
FoxOl antibody and used as a second negative control. 

ChlPed DNA was processed for sequencing by Solexa/ 
Illumina Genome Analyzer II, which generated a data set 
of 5-6 million reads in each group (Figure 1A). The peaks 
were mapped onto the mouse reference genome mm9 in 
the UCSC Genome Browser Database. FoxOl -binding 
peaks were identified by SEAN (Supplementary 
Materials and Methods SI) with the threshold P- value 
1 x 10 -5 . We found 583 peaks in FoxOl -enriched chroma- 
tin compared with IgG in WT (WT-FoxOl/WT-IgG), 
designated FoxOl /IgG, whereas 460 peaks were identified 
in FoxOl -enriched chromatin in WT compared with that 
in 1-FoxOl KO (WT-FoxOl/l-FoxOl KO-FoxOl), 
designated WT/KO (Figure IB). Notably, the overlapping 
peaks were 401, accounting for as much as 69% of peaks 
in FoxOl/IgG and 87% of peaks in WT/KO. The genomic 
locations of the 401 peaks are shown in Supplementary 
Table S9. We used the overlapping peaks as a reference for 
further analysis of genome-wide FoxOl -binding sites. We 
randomly selected three representative examples of 
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Genotype 


Antibody 


Total reads 


WT 


FoxOI 


6,680,589 


WT 


igG 


6,024,711 


KO 


FoxOI 


5,779,448 


401 

182 69%/87% 


59 


Fox01/lgG, N=583 



WT/KO, N=460 




Figure 1. ChlP-seq analysis for FoxOI binding in hepatic chromatin. 
(A) Total ChlP-seq reads are shown in WT or 1-FoxOl KO (KO) 
chromatin enriched with an anti-FoxOl or IgG. (B) Peaks were 
found by SEAN, as described in 'Supplementary Materials and 
Methods SI'. FoxOI -binding peaks in WT FoxOI /WT IgG (FoxOI/ 
IgG) were compared with those of WT FoxOI /KO FoxOI (WT/KO) 
and overlapping peaks were identified. 



ChlP-seq peaks, as shown in Figure 2. The location of the 
peaks for Pdk4 and Pckl-a encompasses the previously 
known FoxOI -binding sites (3,27), whereas the peaks for 
Elk4 and Pckl-b are newly found as putative FoxOI - 
binding sites in the current studies. 

To evaluate the quality of genomic FoxOI -binding 
peaks, we generated density map plots, where each read 
in a peak was displayed as a green spot on a dark back- 
ground horizontally spanning 5kb upstream and down- 
stream from the center of the peak (Figure 3A). 
Therefore, the more reads you have, the brighter the 
spots are. The reads in all of the peaks were presented 
vertically. As shown in Figure 3A (top panel), there was 
a clearly visible vertical bright green line at the center of 
the FoxOI -enriched DNA sample in WT (column 2), 
whereas there was no noticeable line at the center in 
WT-IgG (column 1) and 1-FoxOl KO-FoxOl 
(column 3), both of which were used as negative 
controls. When the peaks in WT-FoxOl were normalized 
against WT-IgG or 1-FoxOl KO-FoxOl, a more distinct- 
ive vertical line was found at the center in columns 4 and 
5, respectively. As a negative control, lOkb of genomic 
regions were randomly selected, and reads in the peaks 
were plotted on a density map (Figure 3A, bottom 
panel). There were no vertical lines visible in the density 
maps generated from random peaks in every case. The 
peak reads in the data set of Figure 3A were plotted in a 
two-dimensional array where the x-axis represents the 
location of reads relative to the center of the peak 
(marked as 0) 5kb upstream and downstream, and the 
j-axis represents the cumulative number of reads in that 
specific location (Figure 3B). Reads from all peaks were 
summed up and plotted. The maximum number of reads 
was set at 100. There was a clear cumulative spike shown 
in WT-FoxOl; however, no visible spike was generated 
from the random peaks of WT-IgG and 1-FoxOl-KO 
FoxOI. The cumulative peak maps that were produced 
from samples normalized against negative controls, 
FoxOI /IgG and WT/KO, showed a higher spike than 
that of WT-FoxOl without normalization. Overall, the 
results in Figure 3 clearly demonstrate that there were 



no distinctive peaks present in two negative controls and 
that the FoxOI peaks are highly specific. 

The analysis of peak distribution demonstrates that the 
preferential binding sites are localized to the intergenic 
regions (37.4%) and introns (34.5%), whereas 15.4% of 
binding sites were mapped within the proximal promoter 
regions. The remaining binding sites were found in 5' 
untranslated regions (8.4%), exons (3.4%) and 3' untrans- 
lated regions (1.0%) (Figure 4A). The analysis of peak 
distribution relative to a known TSS revealed a signifi- 
cant enrichment of FoxOI -binding sites around TSS 
(Figure 4B), whereas the distribution of random peaks 
was spread out relatively even throughout the genome. 
The search for enriched motifs within the FoxOI -binding 
peaks using MEME (21) revealed a sequence element, 5' G 
TAAACA 3', as the highest-scoring motif that was found in 
362 (90%) of 401 peaks with an E = 3.1e-150 (Figure 4C). 
It turns out that this motif closely corresponds to a consen- 
sus site that has been previously known as a FoxOI -binding 
site, GTAAACAA (5,28). There were no other variants of 
the consensus motif found in the MEME analysis. It has 
been shown that FoxOI is associated with other transcrip- 
tion factors to regulate its target gene expression (29). 
To identify enriched binding sites for co-regulatory tran- 
scription factors that are associated with FoxOI, we 
screened the sequence spanning 250 bp upstream and 
downstream from the FoxOI motifs. The motif search 
revealed transcription factor-binding sites for hepatocyte 
nuclear factor 4 (HNF-4), estrogen-related receptor alpha 
(ERRoc), CCAAT/enhancer binding protein beta (CEBPP) 
and glucocorticoid receptor (GR) (Supplementary 
Table S3), in addition to FoxOI. 

FoxOI has been shown to bind to its consensus DNA 
sites with higher affinity when the binding site is flanked 
with its reverse complementary sequence and a specific 
sequence, 5' ACAACA 3' (30). Therefore, we investigated 
whether there are any enriched specific sequences in the 
flanking regions of the FoxOI -binding sites and whether 
there are multiple FoxOI -binding sites present in a peak. 
We found that the majority of peaks contain one FoxOI - 
binding site (81%), whereas 17% of peaks harbor 
two of them. The remaining peaks contain three or 
four FoxOI -binding sites in a peak (Supplementary 
Figure S6A). The number of nucleotides between the 
FoxOI -binding sites ranged from 1 to 182, without any 
specific enrichment in the number of nucleotides between 
them (Supplementary Figure S6B). We did not observe 
any enrichment of particular sequences in the flanking 
regions in our study. 

To evaluate how the genome-wide binding analysis cor- 
relates with global expression profiling of genes regulated 
by FoxOI, we performed GSEA (23). Nearest genes 
within 20 kb of FoxOI -binding peaks were identified in 
the ChlP-seq analysis (Supplementary Table S8). To 
obtain the gene expression profiling for GSEA, we 
carried out genome-wide mRNA expression analysis in 
livers of WT mice either under fed or fasted conditions. 
As FoxOI is mostly active in regulating its target genes 
during fasting, the genes identified by the expression 
profiling assays most likely represent physiologically 
relevant putative FoxOI targets. Additionally, 
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Figure 2. Representative view of a ChlP-seq peak. FoxOl -binding peaks were mapped onto the mouse reference genome mm9 and visualized using 
the UCSC Genome Browser. Arrows indicate peaks associated with genes: (A), Elk4\ (B), Pckl; (C), Pdk4. WT = wild-type; KO = 1-FoxOl KO; 
Fx = FoxOl. 



genome-wide mRNA expression analysis was performed 
in livers of WT and 1-FoxOl KO mice under fasted con- 
dition to prepare another set of FoxOl -regulated gene 
expression profiling. Differentially expressed genes in 
each of the mRNA expression analyses were rank-ordered 
and input into GSEA along with the data set of nearest 
genes in our ChlP-seq analysis. The outcome was an en- 
richment plot, shown in Figure 5 and Supplementary 
Figure S3. The results in Figure 5 demonstrate a signifi- 
cant association of the high density of ChlP-seq genes in 
the running sum plot with the genes that were differen- 
tially expressed with higher fold changes between fed and 
fasted conditions, indicating that the two data sets are 
significantly correlated. The enrichment score (ES) 0.6 
was obtained by running GSEA with 8350 as number of 
permutation with P- value less than 1. 818181 E- 5. A similar 
significant correlation was observed between the two data 
sets from ChlP-seq and from gene expression profiling 
obtained from WT versus KO under fasted conditions 
(Supplementary Figure S3). Fasting triggers adaptive 



changes in metabolic processes primarily for carbohy- 
drates, lipids and amino acids in the liver (31). We show 
in Supplementary Table S4 examples of ChlP-seq genes in 
the metabolic processes, which have been significantly 
correlated with genes in the mRNA profiling. Functional 
categories of genes were obtained by GO analysis, as 
described below, to analyse ChlP-seq genes. Insulin-like 
growth factor-binding protein 1 (Igfbpl) is shown as a 
positive control, as the expression of the gene is directly 
regulated by FoxOl and changed during fasting (32,33). 

To analyse the putative biological function of FoxOl 
ChlP-seq genes, we performed GO analysis using the 
National Institutes of Health Database for Annotation, 
Visualization and Integrated Discovery GO (22). When 
the genes nearest within 20 kb of the FoxOl -binding 
peaks were analysed by broad biological categories, the 
metabolic process (GO:0008152, 54%) and cellular 
process (GO:0009987, 47%) turned out as highest 
enriched terms (Supplementary Table SI). When we per- 
formed the analysis for more specialized biological 
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Figure 3. Density map plots of FoxOl ChlP-seq peaks. (A) The top 
panel represents the read distribution of peaks in ChlP-seq for FoxOl, 
whereas the bottom panel represents the read distribution of random 
peaks. (B) The top and the bottom plots represent cumulative reads of 
FoxOl ChlP-seq peaks and those of random peaks, respectively. The 
maximum number of read is set to 100. 
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Figure 4. Genome-wide distribution of FoxOl -binding sites. 
(A) Genomic distribution of FoxOl -binding sites relative to mouse 
RefSeq genes is shown. The promoter regions are defined as 2kb of 
5' flanking sequence. (B) Peak distribution relative TSS is shown. The 
peak distribution of FoxOl binding is shown in red, whereas the dis- 
tribution of random peaks is shown in blue. (C) The top-scoring motif 
present in the FoxOl -binding peaks is shown by MEME motif analysis. 
E = 3.1e-150. 



Enrichment Plot 




0 5 10 15 20 x10 3 

Rank of genes in fasted/fed microarray 



Figure 5. Gene set enrichment analysis. All genes (n = -2.5 x 10 4 ) m 
microarray expression profiling from fed versus fasted mouse livers 
were ranked by fold difference and P-value and expressed on the 
x-axis. This data set was compared with the gene list of nearest genes 
in = 296) identified in ChlP-seq peaks located within 20 kb of a known 
gene. The plot of running sum for genes in ChlP-seq data is shown at 
the bottom with the leading-edge subset. The location of maximum 
enrichment score indicates a high correlation of the two data sets. 



processes, we found that FoxOl ChlP-seq genes are sig- 
nificantly associated with metabolic processes for 
carboxylic acids, steroids, fatty acids and amino acids 
(Supplementary Table S2). Some genes were also 
categorized for biological processes involved in regulation 
of DNA replication and gene expression/transcription. 
Unexpectedly, with a high fold enrichment score, we also 
observed a category of retinoid (vitamin A) metabolic 
processes where dehydrogenase/reductase SDR family 
member 9 (Dhrs9), retinol saturase (Retsat), Rbpl and 
Rdh8 were included. Interestingly, retinoids have been 
shown to play an important role in hepatic gluconeogenesis 
(13-15). Given that FoxOl plays a critical role in hepatic 
gluconeogenesis, we were prompted to investigate the in- 
volvement of FoxOl in retinoid metabolism further and 
its association with hepatic gluconeogenesis. 

We first confirmed the FoxOl binding to the genomic 
region of the genes in retinoid metabolic processes. Gene- 
specific manual ChIP was performed with an anti-FoxOl 
or IgG. The chromatin peak regions bearing the FoxOl - 
binding motif were amplified by qPCR with gene specific 
primers. Additionally, we randomly picked several 
FoxOl -binding peak regions and the enrichment of 
FoxOl was verified. The peaks for Elk4, Pckl and Pdk4 
are visualized in Figure 2 as examples. The results dem- 
onstrate significant enrichment of FoxOl in the chromatin 
peak regions associated with the genes in retinoid metab- 
olism as well as others selected randomly (Table 1). We 
used Pdk4 and Pckl -a as positive controls for FoxOl 
binding and Rpl32 was used as a negative control. We 
noticed minor peaks for Pckl (Figure 2), which were not 
identified as significant peaks by our peak finding 
program. qPCR confirmation of FoxOl binding to the 
regions indicated no significant enrichment of FoxOl 
binding (data not shown). 

Because FoxOl activity is increased during fasting, 
regulating its target genes, we wondered whether the ex- 
pression of the putative FoxOl target genes in retinoid 
metabolism is altered during fasting. As shown in 
Supplementary Figure S4, there was a significant elevation 
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Table 1. Gene-specific ChIP analysis was performed with an 
anti-FoxOl or IgG 
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Twelve FoxOl -binding regions were randomly picked from ChlP-seq 
data, and the chromatin regions that harbor the FoxOl -binding site 
were amplified by qPCR using gene-specific primers. Results are 
expressed as fold change in comparison of DNA enriched by an 
anti-FoxOl with that of IgG. Data represent the mean of triplicates. 
a,b Indicate the peaks Pckl-a and Pckl-b, shown in Figure 2, 
respectively. 



of Retsat expression by fasting, consistent with a previous 
report (34). The expression of Dhrs9 and Rdh8 was not 
significantly changed. We observed a mild but statistically 
significant reduction of Rbpl expression by fasting. 
Expression of Pckl and sterol regulatory binding protein 
lc, a key regulator of hepatic lipogenesis (35), was shown 
as positive and negative controls for fasting, respectively. 
There was no significant change in FoxOl mRNA levels. 

If FoxOl is an important player in retinoid metabolism, 
deletion of FoxOl should result in impaired retinoid sig- 
naling in regulating its downstream signaling targets. To 
test this idea, we knock-downed FoxOl in primary hep- 
atocytes by an shRNA adenovirus targeting FoxOl 
(FoxOl shRNA) and analysed mRNA expression of 
genes for Pckl and Pdk4 in response to retinol or atRA, 
as Pckl and Pdk4 have been shown regulated by retinoids 
(13,14,16-18). Efficient knockdown of FoxOl by FoxOl 
shRNA is shown in Supplementary Figure S5 along with 
significantly reduced expression of Pckl, a known FoxOl 
target. In line with our ChlP-seq analysis, the expression 
levels of genes for Retsat, Rbpl, Dhrs9 and Rdh8 in 
retinoid metabolism were significantly lower in hepato- 
cytes where FoxOl expression was disrupted. As shown 
in Figure 6, expression of Pckl was significantly induced 
by retinol and atRA, consistent with the presence of 
retinoic acid response elements in the promoter region 
(13,16). PDK4 phosphorylates and inhibits the pyruvate 
dehydrogenase complex, promoting the conversion 
of pyruvate to lactate that is used for hepatic gluco- 
neogenesis. We also observed markedly increased Pdk4 
expression by retinol and atRA. In contrast, the gene 
encoding glucose-6-phosphatase showed no significant 
response to retinol and atRA. Notably, deletion of 
FoxOl resulted in a significantly blunted induction of 
Pckl and Pdk4 in response to retinol and atRA, suggest- 
ing that FoxOl is required for full response of the genes 
to retinoids. 




G6pc 



Ctrl Rol atRA Ctrl Rol atRA Ctrl Rol atRA 



■ UN RNAi □ FoxOl shRNA 

Figure 6. FoxOl is required for full response of Pckl and Pdk4 to 
retinoids. Primary mouse hepatocytes were infected with an adenovirus 
expressing FoxOl shRNA or un-specific RNAi for 24 h. Hepatocytes 
were treated with retinol (20|iM) or atRA (20|iM) for 6h. mRNA 
levels were measured by qPCR and were normalized for ribosomal 
protein Rpl32 mRNA levels. Results are expressed as fold change 
relative to those of control. The mean values obtained from triplicates 
in each group as shown with error bars. *P<0.05. Ctrl = control; 
rol = retinol; atRA = all-trans RA. 



In addition to FoxOl, FoxA2, a member of the FoxA 
subfamily of forkhead box proteins, plays important roles 
in hepatic metabolism (36). FoxA2 has been shown to be 
involved in hepatic glucose homeostasis by regulating 
genes for Pckl, G6pc and tyrosine aminotransferase 
(37-39). FoxA2 also has been shown to be involved in 
hepatic lipid metabolism (40). To investigate whether 
FoxOl and FoxA2 share some aspects of their biological 
functions in the liver or whether they have their own dis- 
tinctive function, we compared our FoxOl ChlP-seq data 
with those of FoxA2 ChlP-seq. Examples of peaks that 
represent a dual target by FoxOl and FoxA2, and a dis- 
tinctive target by FoxOl or FoxA2, are shown in 
Supplementary Figure S7. Tmeml40 is a putative dual 
target by FoxOl and FoxA2, whereas Fig la is a 
putative FoxOl -specific target. Bcas3 is shown as a 
putative Fox A2- specific target. We found that 58.6% of 
FoxOl ChlP-seq genes were co-targeted by FoxA2, 
whereas the remaining FoxOl ChlP-seq genes turned 
out as FoxOl -specific genes (Supplementary Figure S8). 
A GO analysis suggests that FoxOl and FoxA2 share 
their roles in metabolic processes for carboxylic acids, 
lipids and steroids (Supplementary Table S5). 
Interestingly, the category of vitamin A metabolic 
processes appeared as a functional group that is shared 
by FoxOl and FoxA2. In addition to their shared func- 
tions, we found that FoxOl and FoxA2 may have distinct 
function in DNA metabolism and liver development, 
respectively (Supplementary Tables S6 and S7). 



DISCUSSION 

The current study provides high-resolution mapping of 
FoxOl -binding sites throughout the mouse genome in 
hepatic chromatin. We identified 401 FoxOl -binding lo- 
cations on a genome- wide level. To identify functionally 
significant FoxOl -binding locations, we used two inde- 
pendent negative controls in our ChlP-seq study. One 
was an IgG control for non-specific binding of FoxOl 
antibody in chromatin from WT mice, as an IgG control 
is used widespread in conventional ChlP-seq studies. 
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Because there are four members of the mammalian FoxO 
family that share a conserved DNA-binding domain (5), 
we used 1-FoxOl KO mice to establish the specificity of 
FoxOl against other FoxOl proteins. Chromatin from 
1-FoxOl KO mice was immunoprecipitated with an 
anti-FoxOl and used as a second mock control. This 
approach would eliminate the possibility that ChlP-seq 
data could be driven erroneously from the cross-reactivity 
of anti-FoxOl with other FoxO proteins, supporting that 
the peaks found in our studies are likely representative of 
binding sites specifically for FoxOl. Notably, the density 
maps of negative controls, IgG in WT and anti-FoxOl in 
1-FoxOl KO, (Figure 3 A, columns 1 and 3 in the top 
panel) show no visible distinct line at the center that is 
generated from peaks otherwise shown in ChIP with an 
anti-FoxOl in WT (Figure 3 A, columns 2, 4 and 5 in the 
top panel). 

The analysis of FoxOl -binding peak distribution dem- 
onstrates that the majority of peaks were located distant 
from TSS. This trend of genome-wide transcription factor 
binding is consistent with other studies (41,42). In particu- 
lar, in RAW264.7 macrophages stably over-expressing 
FoxOl, only 5.23% of FoxOl peaks are found in the 
proximal regions (43). Similarly, 6% of binding sites for 
FoxA2 in mouse livers were mapped within 1 kb from TSS 
(20). Although intergenic binding sites of some transcrip- 
tion factors have been suggested to mediate biological 
functions of the factors (44,45), it remains largely elusive 
what the precise biological functions of transcription 
factor bindings in non-proximal promoter regulatory 
regions would be. 

We compared our ChlP-seq data with global analysis of 
gene expression regulated by FoxOl in mouse livers 
(Figure 5 and Supplementary Figure S3). The high 
degree of correlation between the two data sets demon- 
strates the functional relevance of our FoxOl ChlP-seq 
genes. During fasting conditions in liver, FoxOl plays a 
significant role in glucose production by regulating 
gluconeogenesis, whereas it functions in suppressing lipid 
accumulation (3,4,46). Consistent with the role of FoxOl 
during fasting, GO analysis of ChlP-seq genes found in 
GSEA showed significant enrichment of functional groups 
of genes in metabolic processes for carboxylic acids, fatty 
acids and amino acids. Supplementary Table S4 shows 
examples of putative novel FoxOl targets in the metabolic 
processes along with known FoxOl targets, Pckl and 
Igfbpl. As the three metabolic processes are shown most 
significantly affected during fasting (31), the data support 
the relevance of FoxOl ChlP-seq genes in hepatic metab- 
olism during fasting. 

In line with these findings, motif search within the FoxOl- 
binding peaks by MEME revealed 5 / GTAAACA3 / that 
closely matches a consensus FoxOl -binding site (28). It is 
enriched in 90% of the peaks with the highest enrichment 
score (E = 3.1e-150) (Figure 4C). FoxOl has been shown to 
cooperate with multiple transcription factors to regulate its 
target gene expression (29). The co-enriched transcription 
factor-binding sites (Supplementary Table S3) found in 
our studies are consistent with previous studies that 
FoxOl is associated with HNF-4, ERRoc, CEBP(3 and 
GR. While FoxOl synergizes with HNF-4 to induce 



transcription of G6pc, it plays a role in modulating the 
activity of ERR and peroxisome proliferator-activated 
receptor gamma coactivator 1 to regulate Pdk4 expression 
(47,48). FoxOl functions cooperatively with CEBPP and 
CEBPoc, an isoform of CEBP(3, to regulate Pckl expression 
in perinatal liver (49). Additionally, FoxOl serves as an ac- 
cessory factor for GR-mediated transcription of Pckl, 
Igfbpl and Pdk4 possibly by recruiting and directly interact- 
ing with p300/CREB-binding protein (CBP), which also 
binds to GR (50-52). These co-regulatory transcription 
factors play important roles in hepatic metabolism, support- 
ing the critical biological function of FoxOl for the meta- 
bolic network in the liver. 

It has been shown that the binding affinity of FoxOl to 
its consensus site is enhanced by the presence of its reverse 
complementary sequence and a specific 6-mer sequence in 
the flanking regions of the binding site. This increased 
binding of FoxOl has been shown to enhance transcrip- 
tional activity of FoxOl in regulating its target gene ex- 
pression (30). Interestingly, the results from our studies 
that the majority of peaks contain one FoxOl -binding 
site suggest that in vivo, FoxOl is able to bind to a 
single DNA-binding site. This is in line with a previous 
study showing that a single consensus motif is sufficient 
for FoxOl binding (28). Among the peaks containing two 
binding sites, we observed that 42% of duplicate FoxOl- 
binding sites consist of a binding site and its reverse com- 
plementary copy in the flaking regions, and that the re- 
maining duplicate FoxOl -binding sites are arranged in a 
direct repeat manner. Therefore, it is likely that the 
multiple DNA-binding sites can be either in a direct- or 
inverted-repeat manner. 

GO analysis of FoxOl and FoxA2 overlapping genes 
demonstrates that FoxOl and FoxA2 may share meta- 
bolic processes for carboxylic acids, lipids and steroids. 
Indeed, FoxOl and FoxA2 play roles in hepatic glucose 
homeostasis by controlling gluconeogenic gene expression 
(3,4,37,39). Additionally, both transcription factors have 
been shown to be involved in lipid metabolism. FoxOl 
controls hepatic lipid accumulation by increasing expres- 
sion of gene for microsomal triglyceride transfer protein to 
promote triglyceride (TG) secretion during fasting (46), 
whereas FoxA2 increases fatty acid oxidation and TG 
secretion, reducing hepatic TG accumulation (40). 
Interestingly, we observed that in GO analysis of 
FoxOl -only or FoxA2-only genes, different sets of genes 
were categorized in some of the same functional groups 
found in the GO analysis of overlapping genes. Therefore, 
it is possible that FoxOl and FoxA2 may share some 
metabolic pathways through regulating their common 
targets as well as their distinct targets in the same 
pathways. Nevertheless, it is likely that they may have 
their distinct functions, shown in Supplementary Figures 
S6 and S7. 

Starvation triggers a complex array of physiologic and 
biochemical pathways to adapt changes in energy metab- 
olism. In liver, glucose production is increased by glyco- 
genosis and gluconeogenesis from non-carbohydrate 
precursors to meet the need for glucose as an energy 
source, especially for the brain and nervous system, 
whereas there is increased fatty acid oxidation and 
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ketone body production. Amino acids serve as non- 
carbohydrate precursors for gluconeogenesis, and there 
is a significantly increased amino acid catabolism during 
fasting. GO analysis of the FoxOl ChlP-seq genes 
demonstrated that the ChlP-seq genes are most highly 
enriched in broad categories of metabolic processes. In 
particular, they were significantly associated with meta- 
bolic processes for carboxylic acids, steroids and fatty 
acids, which is in accordance with known functions of 
FoxOl shown in previous studies. A functional category 
of metabolic processes for cellular amino acid metabolism 
was also observed. Interestingly, there was an additional 
category of gene cluster that is involved in retinoid 
metabolic processes. In the current report, we demonstrate 
the possibility of a previously unappreciated role of 
FoxOl in retinoid metabolic processes and suggest that 
FoxOl may couple retinoid metabolism to hepatic 
gluconeogenesis. 

It is generally thought that retinoids exert their roles in 
biological functions by their active metabolite, atRA, a 
ligand for RARs. The physiological synthesis of atRA is 
tightly controlled through processes involving retinoid- 
binding proteins and retinoid-metabolizing enzymes (12). 
RetSat catalyzes the saturation of all- trans retinol to (R) 
all-trans 13,14-dihydroretinol (dhretinol). Interestingly, 
dhretinol has also been shown to bind to RARs, 
regulating RAR-dependent gene expression (53). It 
remains unknown whether the concentrations of retinoid 
metabolites are being changed during the period of tran- 
sition from the fed to fasted states in liver. Based on our 
study, however, it would be possible that altered expres- 
sion of genes, including Retsat and Rbpl, during fasting 
might influence the availability of retinoid metabolites for 
RARs, which in turn affects their downstream signaling 
targets. Alternately, it is possible that altered retinoid me- 
tabolites could function independently of atRA. Indeed, 
recent studies have shown that parental vitamin A 
(retinol) and retinaldehyde exert their action through sig- 
naling that is not mediated by RARs/retinoid X receptors 
(54,55). It would be interesting to investigate further how 
the alterations of RetSat and RBP1 modulate their down- 
stream signaling pathways. 

Regulation of the Pckl gene by retinoids has been well 
documented in vivo and in vitro. Two retinoic acid 
response elements that mediate retinoic acid responsive- 
ness have been identified on the Pckl promoter (13,16). 
Pckl encodes the rate-limiting enzyme in hepatic 
gluconeogenesis, and the expression of the gene is 
induced by fasting. Importantly, this fasting-mediated in- 
duction of Pckl is blunted by vitamin A deficiency (14). 
Additionally, increased retinol-binding protein-4 (RBP-4), 
a transport protein for retinol in the circulation, causes 
insulin resistance and elevated hepatic Pckl expression 
(56). Similar to Pckl, Pdk4 expression is markedly 
reduced in vivo by vitamin A deficiency (17,18). These 
findings suggest that retinoids play a significant role in 
regulating hepatic glucose metabolism. Because aberrant 
expression of Pckl and Pdk4 is associated with insulin 
resistance in type 2 diabetes, tight control of retinoid 
metabolic processes might be needed to maintain glucose 
homeostasis. 



Our current studies show that the deletion of FoxOl in 
primary hepatocytes resulted in blunted induction of Pckl 
and Pdk4 in response to retinol and atRA, suggesting that 
FoxOl is required for full response of the genes to retin- 
oids. Based on our studies, we propose that during feeding 
conditions where there is increased circulating insulin, 
phosphorylated FoxOl in the cytoplasm would be inef- 
fective in activating its target genes, including those in 
retinoid metabolism, which would result in diminished 
effects of retinoids on its downstream targets. Under 
fasting conditions when the circulating levels of insulin 
go down, nuclear FoxOl would regulate its target genes, 
including genes in retinoid metabolism. This would in 
turn regulate downstream targets of retinoids such as 
Pckl and Pdk4, to supply glucose needed during fasting. 
In support of our notion, retinoid-mediated regulation of 
Pckl gene expression is significantly blocked by insulin 
(15). Hyperglycemia is the hallmark of type 2 diabetes, 
which is highly associated with insulin resistance. Under 
the insulin-resistant status, overactive FoxOl promotes 
hepatic glucose production. We propose that alterations 
in gene expression in retinoid signaling by overactive 
FoxOl might contribute in part to the hyperglycemia 
associated with insulin resistance. There is growing 
evidence that retinoids play a role in energy expenditure 
and adiposity (55,57,58). Consistently, altered retinoid sig- 
naling pathways have been implicated in type 2 diabetes 
and insulin resistance. Future studies will be warranted to 
explore the involvement of FoxOl in retinoid signaling 
and its crosstalk with insulin in diabetes. 
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